arXiv: 1509.00799v2 [hep-lat] 25 Feb 2016 


ADP-15-30/T932 
DESY 15-158 
Edinburgh 2015/20 
Liverpool LTH 1055 
Eebruary 29, 2016 


QED effects in the pseudoscalar meson sector 

R. Horsley®, Y. Nakamura^, H. Perlt®, D. Pleiter^, 

P. E. L. Rakow®, G. Schierholz-^, A. Schiller®, R. Stokes^, 

H. Stiiben^, R. D. Young^ and J. M. ZanottP 

- QCDSF-UKQCD Collaboration - 


“ School of Physics and Astronomy, University of Edinburgh, 
Edinburgh EH9 3ED, UK 

^ RIKEN Advanced Institute for Computational Science, 

Kobe, Hyogo 650-0047, Japan 

Institut fiir Theoretische Physik, Universitat Leipzig, 

04109 Leipzig, Germany 

Jiilich Supercomputer Centre, Eorschungszentrum Jiilich, 

52425 Jiilich, Germany, 

Institut fiir Theoretische Physik, Universitat Regensburg, 

93040 Regensburg, Germany 

® Theoretical Physics Division, Department of Mathematical Sciences, 
University of Liverpool, Liverpool L69 3BX, UK 

Deutsches Elektronen-Synchrotron DESY, 

22603 Hamburg, Germany 

® CSSM, Department of Physics, University of Adelaide, 
Adelaide SA 5005, Australia 

^ Regionales Rechenzentrum, Universitat Hamburg, 

20146 Hamburg, Germany 


February 29, 2016 


Abstract 

In this paper we present results on the pseudoscalar meson masses from a fully 
dynamical simulation of QCDJ-QED, concentrating particularly on violations of 
isospin symmetry. We calculate the splitting and also look at other isospin 

violating mass differences. We have presented results for these isospin splittings 
in [T]. In this paper we give more details of the techniques employed, discussing 
in particular the question of how much of the symmetry violation is due to QCD, 
arising from the different masses of the u and d quarks, and how much is due 
to QED, arising from the different charges of the quarks. This decomposition is 
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not unique, it depends on the renormalisation scheme and scale. We suggest a 
renormalisation scheme in which Dashen’s theorem for neutral mesons holds, so 
that the electromagnetic self-energies of the neutral mesons are zero, and discuss 
how the self-energies change when we transform to a scheme such as MS', in which 
Dashen’s theorem for neutral mesons is violated. 


1 Introduction 

Lattice calculations of the hadronic spectrum are now reaching a precision where it is 
essential to resolve the influence of isospin breaking effects. These have two sources, a 
QCD effect arising from the fact that the u and d quarks have different masses, and 
an electromagnetic effect due to the u and d having different electric charges. The two 
effects are comparable in magnitude, so a reliable calculation of isospin breaking requires 
simulating both the gluon and photon gauge helds. 

Lattice studies of electromagnetic effects in the pions go back to [2]. In recent years 
the interest in QCD-I-QED has grown, and the pace of work accelerated [SHH]. 

We are carrying out simulations in QCD-I-QED [T]. Both gauge theories are fully 
dynamical, so that the electrical charges of sea-quark loops are included via the fermion 
determinants. We use a non-compact action for the photon held. The calculations are 
carried out with three clover-like quarks. Details of the lattice action will be given in 
section lU and can be found in [HIlQ] . 

In the real world, with aEM = 1/137, electromagnetic effects on masses are at the 1% 
level, or smaller. This would make them hard to measure on the lattice. Therefore we 
simulate with a QED coupling stronger than in real world, so that we can see effects easily, 
and then scale back to physical aEM- The simulations are carried out with I3qed = 0.8, 
equivalent to = 1.25, aEM = e^/(47r) 0.10. We will see that this is a good choice, 

electromagnetic signals are clearly visible, much larger than our statistical errors, but we 
are also in a region where they still scale linearly in e^, and we do not need to consider 
higher-order terms. 

We generate conhgurations with dynamical u, d and s quarks, and then increase our 
data range by carrying out partially quenched calculations, with valence u, d, s quarks 
having different masses from the quarks used in the generation of the conhgurations. In 
addition to the u, d, s quarks, we also introduce a hctitious n quark, an extra havour with 
electrical charge zero. The n quark is particularly useful for checking that we are in the 
region where electromagnetic effects are still linearly proportional to e^. 

In this work we present results on the pseudoscalar mesons. Our meson propagators 
are calculated from connected graphs only. Because we have no fermion-line disconnected 
graphs, the uh, dJ, ss and nn states do not mix, so we can measure M‘^{uu), M‘^{dd) and 
M‘^{ss). In the real world, these states do not exist, they mix strongly to form the 7^^,r] 
and r]'. Disconnected graphs are responsible for the large mass of the rj', but will have very 
little effect on the mass of the vr^. In this work we do not consider the r] and r]' further, 
but we will need a mass for the 7r°, with wave-function proportional to {uu — dd) /\/2. We 
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use the relation 


M^o ^ i [M^{uu) + M‘^{dd)\ (1) 

which is a very good approximation, with corrections proportional to the small quantity 
{md—rriuY [H]- This issue does not arise for the flavour non-diagonal mesons, tt’*', 
which have no disconnected contribution. 

In the hrst part of this paper, sections [2] to [TJ we discuss theoretical questions. First we 
describe how our constant singlet mass procedure mm can be applied to QCD-I-QED. 
We derive a mass formula for pseudoscalar mesons in this framework. This is all that is 
needed to calculate physical mass splittings, in particular the 7r+-7r° splitting. It also gives 
us the lattice masses for the u, d, s quarks at the physical point, needed to predict mass 
splittings in the baryons. A particularly delicate number is the mass difference 
(or niu/md mass ratio), which is difficult to extract reliably from a pure QCD simulation, 
and is much better dehned in QCD-I-QED simulations. 

We also want to dissect the meson mass into a QCD part and a QED part, to hnd 
the electromagnetic e parameters, which express the electromagnetic contributions to the 
meson masses [13]. We hnd that there are theoretical subtleties in this separation, leading 
to scheme and scale dependence in the result. 

The total energy-momentum tensor is invariant under renormalisation, and so the 
total mass of any hadron is independent of renormalisation scheme and scale. However 
the individual contributions from quarks, gluons and photons are not invariant, they all 
run as the energy scale increases. This is familiar in pure QCD; as the energy scale of 
Deep Inelastic Scattering rises, the momentum fraction carried by quarks decreases, while 
the momentum fraction carried by gluons increases [13]. The physical picture behind 
this effect is well known [15]. As rises the proton is probed with improved spatial 
resolution. A parton perceived as a single quark in a low-Q^ measurement is resolved into 
multiple partons at higher with most of the new partons being gluons. 

We should expect a similar effect in QCD-I-QED, with improved spatial resolution 
revealing more photons, causing a running of energy from quarks to photons, in parallel 
with the running from quarks to gluons seen in QCD alone. 

In QCD-I-QED, each hadron will be surrounded by a photon cloud. As in pure QED, 
the total energy in the cloud will be ultra-violet divergent. Crudely, we can think of two 
components of the cloud. Firstly, there are short wave-length photons, with wave-lengths 
small compared with a hadron radius. These can be associated with particular quarks. If 
we look at the hadron with some hnite resolution the photons with wavelengths shorter 
than this resolution are incorporated into the quark masses as self energies. Secondly, there 
will be longer wave-lengths photons, which can’t be associated with particular quarks. 
These photons must be thought of as the photon cloud of the hadron as a whole, these 
are the photons that we include when we talk of the electromagnetic contribution to 
the hadron mass. We expect to see many more really long wave-length photons (large 
compared to the hadron radius) around a charged hadron than around a neutral hadron. 

Clearly, in this picture, the value we get for the electromagnetic contribution to the 
hadron energy is going to depend on our resolution, i.e. on the scheme and scale that we 
use for renormalising QED. 
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In the final part, section [HI we summarise our lattice results for the 7r+-7r° splitting 
and for the scheme-dependent e parameters, which parameterise the electromagnetic part 
of the meson masses. 

We have already published an investigation into the QCD isospin breaking arising 
from rrid — rUu alone in ini, and the first results of our QCD-I-QED program in [1], which 
we discuss at greater length here. 


2 Extrapolation Strategy 

In pure QCD we found that there are significant advantages in expanding about a sym¬ 
metric point with = rUg = m [TTIIT^ . In particular, this approach simplifies 

the extrapolation to the physical point, and it decreases the errors due to partial quench¬ 
ing. We want to follow a similar approach with QED added, even though the symmetry 
group is smaller (the u quark is always different from the other two flavours because of 
its different charge). 

First we find a symmetric point, with all three quark masses equal, chosen so that the 
average quark mass, 

m = i (m„ + rnd + mQ , (2) 

has its physical value. To do this, we have defined our symmetric point in terms of the 
masses of neutral pseudoscalar mesons 

M‘^{uu) = M\dd) = M^{ss) = M^{nn) = . (3) 

Here is an average pseudoscalar mass, defined by 

A'J = i f + (M: f] (4) 

where * denotes the real-world physical value of a mass. The n is a fictitious electrically 
neutral quark flavour. We have not included disconnected diagrams, so the different 
neutral mesons of ([3]) do not mix. 

We also define the critical for each flavour as the place where the corresponding 
neutral meson is massless 0 

= 0 rUg = 0 . (5) 

Chiral symmetry can be used to argue that neutral mesons are better than charged ones 
for defining the massless point [T6] . 

We then make a Taylor expansion about this point, using the distance from rn as our 
parameter to specify the bare quark masses 

adrria = airria — rn) = -, 

adjjjq = Cl(/iq m) 1 

^The critical k defined in eq. m is the critical k in the -I- md + rag = const surface, i.e. if m„ = 0, 
we must have rud + nig = 3m. The for the chiral point with all three quarks massless will be different. 


( 6 ) 

(7) 
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where niq denotes the simulation quark mass (or sea quark mass), while /ig represents the 
masses of partially quenched valence quarks. Note that keeping the average quark mass 
constant, ([2]), implies the constraint 

6mu + Srrid + SrUs = 0 . ( 8 ) 


In [H] we wrote down the allowed expansion terms for pure QCD, taking flavour 
blindness into account. QCD+QED works very much like pure QCD. Since the charge 
matrix Q is a traceless 3x3 matrix. 


Q = 




(9) 


electric charge is an octet, so we can build up polynomials in both charge and mass 
splitting in a way completely analogous to the pure QCD case. The main difference is 
that we can only have even powers of the charge, so the leading QED terms are ~ e^, 
while the leading QCD terms are ~ 5m. 


One very important point to note is that even when all three quarks have the same 
mass, we do not have full SU(3) symmetry. The different electric charge of the u quark 
means that it is always distinguishable from the d and s quarks. 


3 Meson mass formula 

From these considerations we hnd the following expansion for the mass-squared of an ah 
meson, incorporating both the QCD and electromagnetic terms 

M‘^{ab) = -|- a{6jaa + Sjai,) + c{Smu -l- Sm^ + 5ms) (10) 

+ /3o|(<^^u + + l^2{5^a — 5fib)‘^ 

+ + + el) + + el) + /3f^(e, - e^f 

+ + el5md el5ms) + '^^^{el5fia + el5fib) 

+ - eQ^(5/ia + 5/iQ + li^{el - el)i5f,a - 5/iQ 

+ 7f^{el + el + el){5f^a + 5f^b) 

+ 7f^(ea + eb){eu5mu + ed5md + . 

As well as the terms needed in the constant fn surface we have also included the term 
c{5mu+5md+5ms), the leading term describing displacement from the constant m surface. 
Including this term will be useful when we come to discuss renormalisation and scheme 
dependence, it could also be used to make minor adjustments in tuning. 

The QCD terms have been derived in |TT]. In particular, we discussed the effect of 
chiral logarithms in section V.C. of that paper. Briefly, since we are expanding about a 
point some distance away from all chiral singularities the chiral logarithms do not spoil 
the expansion, but they do determine the behaviour of the series for large powers of 5mq, 
(see for example equation (78) of [TT|). 

We will now discuss briefly the origins of the electromagnetic terms. 
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3.1 Leading order terms 

In what follows we use the following notation: 

= ^/f^QED , Cq = QqC (11) 

where 

Qu = + |) Qd = Qs = — \ ■ ( 12 ) 

The leading order EM terms were written down in [TO] , 

Mluiab) = + el + el) + {el + el) + - e,)^ . (13) 

Upon examination of each of these terms in more detail, we observe that since all of 
our simulations have the same choice of sea quark charges, then even if we vary the sea 
quark masses, {el + el + el) is a constant, and we can simply absorb this term into 
of flTOl) . Hence, the I3q^ term just stands for the fact that measured in QCD+QED 
might be different from measured in pure QCD. As we have tuned our expansion 
point so that the pseudoscalars have the same symmetric-point mass as in pure QCD, the 
/9(f^ for the pseudoscalar mesons will be zero, but we will still have to allow for other 
particles to be different in QCD-I-QED than in pure QCD. 

Now consider fITU]) at the symmetric point, for the case of a flavour-diagonal meson, 
aa. At the symmetric point, nearly all terms vanish because Srrig and 6fiq are zero. In 
addition, the electromagnetic terms simplify because = e^. All we are left with is 

M\aa) =M^+ + el + el) + 2/3^6^ (14) 

at the symmetric point. However, since we have defined our symmetric point by (I3|) . 
equation (ITT)) must give the same answer whether Ca = — ^e, 0 or +|e, so must be 
zero (because it would split the masses of the different mesons, according to the charge 
of their valence quarks). However, having = 0 for the pseudoscalar mesons does not 
mean that this term will also vanish for other mesons, for example the vector mesons. If 
we tune our masses so that the pseudoscalar uu, dd and ss all have the same mass, we 
would still expect to hnd that the vector uu meson would have a different mass from the 
vector dd and ss, because there is no symmetry in QCD-I-QED which can relate the u to 
the other two flavours. 

Finally, we observe that the contribution from is zero for neutral mesons, Ca = et- 
However, this is the leading term contributing to the mass splitting, so it is of 

considerable physical interest. 

3.2 Next Order 

Going beyond leading order, the following higher order terms of the form e’^drriq, e^5fiq 
are possible: 
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• Sea charge times sea mass, 7 ,f^^ 

After imposing the constraints that m is kept constant and e^l + = 0, there 

is only one completely symmetric sea-sea polynomial left, 

e^drriu + e^drud + e^Snis . (15) 

• Valence charge times sea mass 

At this order all polynomials of this type are killed by the m = const constraint. 

• Valence charge times valence mass, 7^^, 7 ®^, 7®^ 

In this case there are three independent allowed terms. One convenient basis for 
the valence-valence terms is 

V ) (Ca Cb) “1“ ^y (('a ^f^b') y (^6) 

though other choices are possible. 

• Sea charge times valence mass, 

The only polynomial of this type is 

i^l + ^d + + Sfib) ■ (17) 

Since (e^ + held constant, this term can simply be absorbed into the 

parameter a of flTOl) . 

• Mixed charge times sea mass, 

At the symmetric point we can not have mixed charge terms (valence charge times 
sea charge), because such terms would be proportional to (e„ + Cd + e^) which is 
zero. However, away from the symmetric point 

(Ca + eb) {cuSniu + CdSnid + J (18) 

is allowed. 

We illustrate the different physical origins of these terms by drawing examples of the 
Feynman diagrams contributing to each of the electromagnetic coefficients in flTOl) . Fig. [H 


4 Lattice setup 

We are using the action 

S = Sg + Sa + S^ + S'^p + S^j,. (19) 

Here Sq is the tree-level Symanzik improved SU(3) gauge action, and Sa is the noncom¬ 
pact U(l) gauge action of the photon, 

Sa = ^I3qed ^ + h) - "4^(^ + z>) - A^(x)f . (20) 
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Figure 1: Examples of the Feynman diagrams contributing to each of the electromagnetic 
coefficients in the meson mass formula (fTOj) . All the graphs have a single photon (wavy 
line), and are all of 0{e^) in the electromagnetic coupling. However, some terms require 
multiple gluons (curly lines), and so have higher order in the strong coupling 


The fermion action for flavour q is 

X ^ H 


-q{x){'y^ + ^^Ul{x - fi)q{x - fi) 

+ 7^li.x)q{x) - \csw^q{.x)a^^F^^{x)q{x) 




( 21 ) 


where is a singly iterated stout link. We use the clover coefficient csw with the value 
computed non-perturbatively in pure QCD, [18]. We do not include a clover term for the 
electromagnetic held. We simulate this action using the Rational Hybrid Monte Carlo 
(RHMC) algorithm [T9] . 

One issue that arises in the simulation of QED is the treatment of constant electromag¬ 
netic background helds. In simulations where the electromagnetic held does not couple to 
the quark determinant these are electromagnetic zero modes, and so need to be handled 
with particular care. In this simulation the sea quarks are coupled to the electromagnetic 
held, and so the action does depend on the background held. However we do still need to 
give special treatment to these modes. We handle constant background helds by adding 
or subtracting multiples of 67r/(eL^) until the background held is in the range 


— Svr < eB^L^ < Svr 
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( 22 ) 











This is the mildest way to keep the background fields under control [20]. This procedure 
leaves fermion determinants unchanged for particles with charges a multiple of e/3. It also 
leaves Polyakov loops unchanged (again, for charges in units of e/3). We are investigating 
the evolution of these background fields in our simulations, and considering what effect 
they have on finite size effects. We plan to report on these studies in a future paper. 

We have carried out simulations on three lattice volumes, 24^ x 48, 32^ x 64 and 48^ x 96. 
The 24^ X 48 calculations show clear signs of finite size effects. The differences between 
32^ X 64 and 48^ x 96 are quite small, leading us to believe that finite size effects on our 
largest volume are under control. In this paper we present results from the two largest 
volumes, which usually are in close agreement. In the few cases where there is a difference, 
we would favour the results from the largest volume, 48^ x 96. 


5 Critical tv 

After several tuning runs we have been carrying out our main simulations at the point 

Pqcd = 5.50 , Pqed = 0.8 , (23) 

Ku = 0.124362 , Kd = Ks = 0.121713 

which lies very close to the ideal symmetric point defined in ([3]) (but with a much stronger 
QED coupling than the real world, aqED = 0.099472 • • •, instead of the true value 1/137). 
At this point the 6mq from the sea quark masses are all zero, but we can still learn abut 
the meson masses by varying the partially quenched valence quark masses, 6fig. 

The flavour dependence of the meson masses is more complicated in QCD+QED than 
in pure QCD. We illustrate some of these differences in the sketch Fig. [2l showing the 
way that the flavour-diagonal mesons depend on the quark mass. As well as the physical 
charge -|-| and — ^ quarks, we also have a fictional charge 0 quark. In QCD-I-QED we 
still have the relationship M^{qq) oc rrig for flavour-diagonal (neutral) mesons, but the 
gradients of the uu, dd, nfi mesons differ. So, in contrast to pure QCD, equal meson mass 
at the symmetric point no longer means equal bare quark mass. The bare mass at the 
symmetric point depends on the quark charge. This situation is illustrated in the left 
panel of Fig. |2l (though the differences between the flavours has been exaggerated for 
clarity). 

We rescale (renormalise) the quark masses to remove this effect, making the renor¬ 
malised quark masses at the symmetric point equal. The situation after renormalising 
in this way is illustrated in the right panel of Fig. [2J All the flavour-diagonal mesons, 
nn, dd, ss and uu now line up, depending in the same way on the new mass which 
we call the “Dashen scheme” mass, for reasons which should become clear later [§. We 
will see that using this quark mass also simplifies the behaviour of the mixed flavour 
mesons, and helps us understand the splitting of a hadron mass into a QCD part and an 
electromagnetic part. 

^Here, to introduce the idea, we just make a simple multiplicative renormalisation. In fact, the mass 
renormalisation matrix is not diagonal, there are also terms which mix flavours. We will include these 
additional terms in section 
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Figure 2: Sketch illustrating the transformation from bare masses (left panel) to Dashen 
scheme masses (right panel). In the left panel all the flavour diagonal mesons have the 
same mass at the symmetric point {Sfig = 0), but have different critical points (Mpg = 0). 
In the Dashen scheme (right panel) we rescale the masses horizontally, so that all the 
critical points are the same. The different mesons now all depend on in the same 
way. 


One way to interpret the behaviour in Fig. [2] is to consider a u and d quark with the 
same bare lattice mass. Since the magnitude of the charge of the u quark is twice as large 
as that of the d quark, it will acquire a larger self-energy due to the surrounding photon 
cloud and hence it will be physically more massive, which is why the mass of the uu 
meson rises more steeply than the dd meson, when plotted against bare mass. By instead 
plotting against the Dashen mass, we have effectively added the extra mass of the photon 
cloud to the quark mass. Two quarks with the same Dashen mass are physically similar 
in mass, and so they form mesons of the same mass, as seen in the right-hand panel of 

Fig. H 

Applying these ideas to our simulations, in Fig. [3] we show how the symmetric 
and critical are determined, using the dd meson as an example, is dehned from the 
point where the partially-quenched meson mass extrapolates to zero, (jS]), while is 
dehned by the point where the ht line crosses Mpg = X^, ([3]). 

We repeat this procedure for the u and n quarks and plot the resulting 1/k^ and 
l/i^sym g^g g function of the square of the quark charges, Q^, in Fig. 0] Here we 

clearly see that in both cases 1/k depends linearly on Q"^. 

Despite appearances, the two lines are not quite parallel. In Fig. [5] we plot the bare 
mass at the symmetric point. 


am 


sym 

Q 


1 

2 ^ 



(24) 


Kg for each Favour is dehned as the point at which the havour-diagonal qq meson becomes 
massless. We see that our data show the behaviour shown in the left-hand panel of Fig. [2l 
with each meson reaching the axis at a diherent point. 

The factors needed to bring the charged bare masses into agreement with the neutral 
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Figure 3: Determination of k'^ and for the d quark, is defined from the point 
where the dd meson mass extrapolates to zero, ([S]), while is defined by the point 
where the £t line crosses Mpg = ([3]). 



Figure 4: 1/k^ (red squares) and (blue circles) plotted against quark charge 

squared, 


bare mass, as in the right-hand panel of Fig. [2l are 

= 1.023, = 1.096 . 

''I'd ''I'S ’ iit'U 


( 25 ) 
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Figure 5: The bare mass at the symmetric point, arriq^'^, as a function of quark charge. 
We see that the bare mass is not constant, there is about a 10% difference between the 
neutral n quark and the u quark. The open red circles show the quark masses after 
renormalising to remove this charge dependence. 


As seen in Fig. |5]this Z factor depends linearly on the quark charge squared. Hence, we 
can write 

= (1 + Kel)5fXq = (1 + KQle^)5iiq , (26) 

for some constant K. By construction, this simplihes the neutral mesons as they will all 
lie on the same line, see Fig. |2l 

In order to investigate the effect on charged mesons, we hrst consider the uu, dd and 
ud (7r+) meson masses plotted as a function of bare quark mass. Fig. [6l We see that in 
this plot the two neutral mesons, uu and dd, lie on different lines. We also observe that 
the 7r+ data do not he on a smooth curve. This is not due to statistical errors (which 
are much too small to see in this plot). It is because the tt’*' meson mass depends both 
on StiIu + Snid, as in pure QCD, but also has a signihcant dependence on StiIu — Snid, 
which causes those mesons containing quarks with very unequal masses to deviate from 
the trend. 

When we now switch to using the Dashen-scheme quark masses in Fig. [7] we see that 
the graph looks signihcantly different. The uu and dd mesons now he on the same straight 
line (this is essentially by construction, since equal Dashen-scheme quark mass equal 
neutral meson mass). More interesting is the fact that the “jiggles” in the mass 
are largely removed by plotting against Dashen-scheme mass, making it much easier to 
estimate the EM shift in the mass. 
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0.06 



- 0.01 0 0.01 0.02 0.03 0.04 0.05 

a((5/i,j^+(5/x^) 


Figure 6: Pseudoscalar Mp^ plotted against bare mass for the (red), uu (blue) and 
dd (black) mesons. The lines simply connect the points. Error bars are small compared 
with the points. Data are from a 32^ x 64 lattice. 



- 0.01 0 0.01 0.02 0.03 0.04 0.05 

a((5/u'^+(5/x°) 

Figure 7: The same data as in Fig.[6l but this time plotted against Dashen-scheme quark 
mass. 
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6 Dashen scheme quark mass formula 

In order to derive an expression for the meson masses in the Dashen-scheme, we start 
with ffTOl) and proceed by absorbing the QED terms for the neutral pseudoscalar mesons 
into the quark self-energy by making the definition 

(5/r^ = 5/ig I 6md + 6ms) + e^dmd + J (27) 

+ + el + el)6iiq + 7f^eg(e„(5m„ ed6md + . 

At present we are neglecting 7 ®^ and because we are working on a symmetric 
background, 6mq = 0, and absorbing into the coefficient a because we only have 
data at one value of (3qed- This means that only the yf^ term is used in calculating 
6fi^, giving a simple multiplicative transformation from bare mass to Dashen scheme 
mass. Most of the other terms in (l27)l represent off-diagonal terms in the quark mass Z 
matrix. There are many more mixing terms possible in QCD-I-QED than in pure QCD, 
but most of them first occur in diagrams with a large number of gluon and quark loops, 
as can be seen in Fig. [H so they are probably rather small. 

Substituting fl7r|) into flTUD we are left with the simpler formula 

M^{ab) = + 6fi^) +/3o^{6m'l +6ml +dml) (28) 

+ + /92(5/if - - e ,)2 

+ yf^(ea - e,)^(<5/if + 5/rf ) + y3^^(e^ - e^)(<5/rf - 5/rf) . 

In fl28ll all the EM terms vanish for neutral mesons (ca = Cb), leaving 

MLutiab) = + a{6^ij^+ 6^1^)+f^oliSml + Sml + Sml) (29) 

+ A )^ + {6f^b)^) + A {6f^a - , 

which clearly has no references to any EM coefficient, or to any charges Cq. Hence, by 
construction, the mass of the neutral pseudoscalar mesons comes purely from the quark 
masses, and has no electromagnetic contribution. The formula simplifies even further if 
we consider a flavour-diagonal meson 

M^(aa) = + 2a6iJ,^ + /3o^{6ml + + 6ml) + 2A(^hf )^ • (30) 

This agrees with what we see in Figs. |2]and[71 with the different flavour-diagonal mesons 
all lying on the same curve when plotted against the Dashen quark mass. 

In the Dashen scheme the electromagnetic contribution to the meson mass is 

M^{ab) = /3f*^(ea-efe)^-t-yf^(ea-efe)^(5pf-Fh/if) (31) 

+ ^i^(el-el)(6^^-6^^), 

while the QCD contribution is 

Mgcniab) = + a{6fi^+ 6fijl’)+l3oli6ml + 6ml + 6ml) (32) 

+ A((5/^f )^ + {6f^m + A(5/^f - 5/if )^ • 
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Dashen’s theorems [21] state that in the limit of an exact SU(3) chiral symmetry, 
the neutral mesons have zero electromagnetic self energy; and that the charged mesons 
electromagnetic self-energies are given by a single constant. Our formulation is such as 
to maintain the vanishing electromagnetic self-energy of the neutral mesons away from 
the chiral limit. The /5|'^ term of our expansion is the generalisation of Dashens result, 
where, in the absence of any strong SU(3) breaking, the electromagnetic self-energy is 
proportional to the charge-square of the meson. The terms involving 7 ®^ therefore encode 
the deviations associated with leading-order SU(3) breaking of the strong interaction, as 
anticipated by Dashen. 


7 Scheme dependence 

We can calculate electromagnetic contributions to the meson masses from fl3Tl) in our 
scheme, but in order to compare our results with those obtained by other groups, we need 
to be able to quote the QED contribution in other schemes, in particular MS. 

To illustrate the issue of scheme dependence, consider the splitting between the 
and mesons. In the real world the K^-K^ splitting comes partly from QED effects, 
and partly from the rrid, rriu mass difference, which we consider to be the QCD part of the 
splitting. The ordering of the physical states, with the heavier than the suggests 
that the quark mass effect dominates, but we expect that there is still a QED contribution 
of comparable magnitude. 

Naively, one might think that this QED contribution may be easily determined by 
performing a simulation with = m^. In this case, there will be no splitting from QCD, 
so the result will give the splitting due to QED alone. In pure QCD, setting 
is unproblematic as equal bare mass implies equal renormalised mass, regardless of scale 
or scheme. However in QED-I-QCD, mass ratios between quarks of different charges are 
not invariant. The anomalous dimension of the quark mass now depends on the quark 
charge; at one-loop 

7 m = + QQy + • • ■ (33) 

so the u mass runs faster than d mass. If = md in one scheme, this will not be true in 
another. This also implies that there is no good way to compare masses at the physical 
with pure QCD masses at = 0. 


7.1 Changing Scheme 


To calculate the electromagnetic part of the meson mass we take the difference between 
the mass calculated in the full theory, QCD-I-QED, and both non-zero) and subtract 
the mass calculated in pure QCD, (e^ = 0): 




(34) 


where e* is the physical value of the electromagnetic coupling, corresponding to aEM = 
1/137. In the full theory the physical quark masses are well dehned: we can £x the three 
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physical quark masses by using three physical particle masses (the 7r°, and would 
be a suitable choice). In the full theory we should use the physical quark masses, m*, but 
we also have to specify which quark masses we are going to use in the pure QCD case, 
(which is, after all, an unphysical theory). Different ways of choosing the will give 

different values for the electromagnetic part of the meson mass. 


One prescription for choosing the quark masses in the (unphysical) pure QCD case is 
to use the neutral meson masses. We could tune by requiring 


el, rrC nij, ml) = 0, 171®°-°, m?®®) 


(35) 


Since the QCD+QED mass matches the QCD mass, this scheme has zero EM contribution 
to neutral pseudoscalars by dehnition. This is our Dashen scheme, discussed above. In 
this scheme, is zero for neutral pseudoscalar mesons, and is given by the simple 
formula flM]) for charged mesons. 

A more conventional choice is to choose m* and mP^^ the same in MS at some 
particular scale. In this case, we are now presented with the task of determining the 
quark masses in a certain scheme (e.g. the Dashen scheme) given hxed MS masses. 
Hence we need to calculate the Dashen quark masses by renormalising from MS to the 
Dashen scheme: 


= Z^{g\el , (36) 


However, since the renormalisation factor Zm depends on both g^ and e^, the Dashen 
mass in pure QCD would not be the same as the Dashen mass in the physical QCD+QED 
theory: 


^QCD = 0, +Q = el) = Yr^{g^, el, ^P)m^{g‘^, el) . (37) 


Zm{g‘^,el,g?) 


Hence the Dashen mass is rescaled by a renormalisation constant ratio which we denote 
Y 

^ rri' 


Now, we know in principle what the QCD mass we should subtract is, it is the mass 
we get by substituting = 0,m^ = into our £t formula. So now it is a mat¬ 

ter of determining the ratio TA in fl57)) To proceed, we note that we already know the 
renormalisation factor from bare lattice mass to Dashen mass, equation (1261) and (|27)) : 


SMlatt^D 
^ m 


EM 

1 + 

a ^ 

1 + cemQI 2.20(9) . 


(38) 


We also need the renormalisation factor from bare lattice mass to MS, which can be 
estimated from lattice perturbation theory [22]. Fortunately, all pure QCD diagrams 
with only gluons and quarks cancel because we are looking at a ratio of Z factors, so the 
leading contribution comes from the 1-loop photon diagram, giving 

_ (p‘(M 

yiatt^MS ^ l + _l|(-61na/r + 12.95241) 

Ibvr^ 

= l + a^MQ? 1.208 . (39) 
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The numerical value in the second line is obtained for /i = 2 GeV and the value of the 
lattice spacing in our simulations, a~^ = 2.9 GeV (see Table [2]). However, the one-loop 
result is not the full answer, there will be higher order diagrams, with one photon plus 
any number of gluons, giving contributions ~ e'^g ^,... To account for these unknown 

terms we add an error ~ ±30% to the coefficient, giving 

= 1 ± aEMQl 1.2(4) . (40) 

Gombining this with fl38|) gives us the conversion factor from the Dashen scheme to MS 
at /i = 2 GeV for our conhgurations (a“^ = 2.9 GeV) 

YD^ms ^ ^ ^ ^ ^ aEMQlr^^'^ . (41) 


We are now ready to write the transformation formula from Dashen scheme to 

(42) 


in MS'. In the Dashen scheme 

r2lD 


[M^r = KJ°. K]°. Kl°) - K|°, K|°, |m 


ii") 


with the same Dashen-scheme quark masses in both terms. In MS 

[M^]^^ = M^ig^^e"^, [m*]^) - M2(/,0, [m, 

where [fhg]^ is given by flTT)) 


lD\ 


[m,f = (1 + aEMQlr^^^^) [m*]^ 


I 


(43) 


(44) 


Taking the difference between fl43D and (1421) gives 


[M=] [M=] ° = MV, 0. K]°. K]°. Kl°) - MV-, 0, K]°, K|°, |mj°) (45) 


which holds for the electromagnetic contribution to any hadron. If we are specifically 
interested in pseudoscalar mesons, we can use the leading order mass formula M‘^{ab) = 
otirua + mb) to give 


= [M^{ab)f- [QlM^aa) + QlM^bb)] . (46) 

This is a rather simple formula, the only difficulty is that at present we only have a rather 
rough value for the constant T. 


8 Lattice Results 

The hrst question to consider is how close our simulation is to the symmetric line, where 
M{uu) = M{dd) = M{ss). We hnd that at the simulation point, M{uu) is about 6% 
heavier than the other two mesons, so we are not quite at the desired point. In Table [T] 
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flavour 

32^ X 64 

48^ X 96 

simulation 

n 

0.1208142(14) 

0.1208135(9) 


d, s 

0.1217026(5) 

0.1217032(3) 

0.121713 

u 

0.1243838(10) 

0.1243824(6) 

0.124362 


Table 1: The k values of the symmetric point, determined from hts to the pseudoscalar 
meson data. 


we show the values determined on our two large-volume ensembles. In our hts we 
make a Taylor expansion about the symmetric point of Table [H not about our simulation 
point. (The displacement is rather small, the difference is in the hfth signihcant hgure.) 

The next question is whether we have the value of rn correctly matched to the physical 
value. This is checked by comparing the averaged pseudoscalar mass squared, (jH), 
with the corresponding baryon scale 

xi, = 1 + (MJ)2 + ■ (47) 

We hnd X^/X^ = 2.79(3), very close to the correct physical value, 2.81, showing that 
our tuning has found the correct fn value very successfully. 

8.1 The splitting of the 7r+ and masses. 

The hrst quantity we wish to consider is the mass difference between the 7r~^ and 7r° 
mesons. Since in this case we are calculating a physically observable mass difference there 
is no scheme dependence in the result. 

First we need to hnd the k values corresponding to the physical quark masses. Since 
we have three quark masses to determine we need three pieces of physical input, we choose 
the masses of the 7r° and the two kaons 

M^o = 134.977 MeV, 

Mko = 497.614 MeV, (48) 

Mk+ = 493.677 MeV 

at aEM = 1/137. This determines the physical point given in Table O We see very close 
agreement between the lattice scale determined on the two lattice volumes. 

Using these quark masses we now have a prediction for the one remaining meson mass, 
the 7r+. Our values on the two lattice spacings are given in Table |3l 

8.2 The e parameters 

The vr^-vr*^ mass splitting that we presented in the previous section is a physically mea¬ 
surable quantity, so it is independent of renormalisation. However, if we now attempt to 
divide our hadron masses into a QCD part and a QED part, as explained earlier, this is 
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32^ X 64 48^ X 96 

a6ml 

a5rnl 

a6m* 

a-VGeV 

-0.00834(8) -0.00791(4) 

-0.00776(7) -0.00740(4) 

0.01610(15) 0.01531(8) 

2.89(5) 2.91(3) 


Table 2: Bare quark mass parameters at the physical point, and inverse lattice spacing, 
dehned from These masses have been tuned to reproduce the real-world 7r°, and 
when aEM = 1/137. 



32^ X 64 48^ X 96 Real World 

M^+ 

A7n-+ — M^J-O 

140.3(5) 139.6(2) 139.570 

5.3(5) 4.6(2) 4.594 


Table 3: The predicted value of the mass, and vr^-vr*^ splitting, in MeV. 


a scheme-dependent concept. When we look with greater resolution we see more short 
wavelength photons, which had previously been counted as part of the quark mass, and 
therefore part of the QCD contribution to the mass. 

The traditional way of expressing the electromagnetic contributions is through the e 
parameters, which measure in units of 

= M/+ - Ml, , (49) 


a natural choice because it is a quantity of a similar origin, and similar order of magnitude. 


The e parameters are dehned by [T3] 

Ml{7^^) = Ml,{g\e^)-Ml,{g\Q) 
MliK^) = MUg^e^) - 
Mlin^) = Ml4g^,e^) - Ml^g^O) 


CttoAt,- , 

[1 -|- — em]^w , (50) 

= [1 + e d- e^o — em]^iT ■ 


eK+ is dehned in this way so that the electromagnetic contribution to the following quan¬ 
tity has a simple expression 


[M|+ - Ml, - Ml+ + Mil = • (51) 

From now on we will neglect the small quantity e^, the QCD contribution to the 
splitting, which comes largely from annihilation diagrams. This is a reasonable assumption 
here since we note that phenomenological estimates for the this QCD contribution are of 
order 0.1 MeV (or 2%) [23], which is within the precision of our present calculation. 

In the Dashen scheme the e parameters are simply, 

efo=0, e^o = 0, ef+ = l, (52) 
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(53) 


with the only non-trivial quantity, e'^, given by 


e 


D 


M2(7r+) 


4 + 


On our two ensembles we hnd 


= 0.38(10) 32^ x 64 , 

= 0.49(5) 48^ x 96 , 


(54) 


which agree within errors. In what follows, we use the 48^ x 96 value in our calculations. 
Using fl46|l to transform these numbers into MS with the scale fi = 2 GeV, we hnd: 

e^o = + ^M^dd)] /A^ = 0.03 ± 0.02 , 

e,+ = I [lM\uu) + lM\dd)]/A^ = 1.03 ±0.02, 

[^M^idd) + iM^ss)] /A, = 0.2 ± 0.1 , (55) 

eK+ = e^+-aEMr^^'^^^[lM^{uu) + lM^{ss)]/A^ = 1.7±0.1, 
e = [|M2(mm) - |M2(dd)]/A^ = 0.50 ±0.06 . 

In all cases we are resolving more photons in MS, and so converting some fraction of the 
quark mass into electromagnetic energy. This has very little effect in the pions because 
both quarks are very light, but a much larger effect in the kaons because the strange 
quark is heavier, and the photon cloud has a mass proportional to the quark mass. 


9 Conclusions 

We have investigated isospin breaking in the pseudoscalar meson sector from lattice cal¬ 
culations of QCD-I-QED. This allows us to look simultaneously at both sources of isospin 
breaking, the quark mass differences, and the electromagnetic interaction, which are of 
comparable importance. 

The physical mass differences between the different particles are directly observable, 
and so must be independent of the renormalisation scheme and scale used. When we try 
to go beyond this, to say what fraction of a hadron’s mass-squared comes from QCD, 
and from QED, this no longer holds — changing our resolution changes the fraction. We 
understand this effect, both formally, in terms of the dependence of the mass renormal¬ 
isation constant on the electromagnetic coupling, and physically, in terms of the quark 
mass gaining a contribution from its associated photon cloud. 

With this understanding, we calculate the electromagnetic contributions to hadron 
masses in the Dashen scheme, which is easy to implement on the lattice, and then convert 
these values into the more conventional MS scheme. 

We are also investigating the isospin violating mass splittings in the baryon sector [1], 
as well as the decomposition of these mass differences into QCD and QED parts, both in 
the Dashen scheme, and in MS. 
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